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Calculating the microscopic dissociation rate of a bound state, such as a classical diatomic molecule, 
has been difficult so far. The problem was that standard theories require an energy barrier over 
which the bound particle (or state) escapes into the preferred low-energy state. This is not the 
case when the long-range repulsion responsible for the barrier is either absent or screened (as in 
Cooper pairs, ionized plasma, or biomolecular complexes). We solve this classical problem by 
accounting for entropic memory at the microscopic level. The theory predicts dissociation rates for 
arbitrary potentials and is successfully tested on the example of plasma, where it yields an estimate 
of ionization in the core of Sun in excellent agreement with experiments. In biology, the new theory 
accounts for crowding in receptor-ligand kinetics and protein aggregation. 
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The rate of escape of a classical particle over an energy 
barrier is a well-posed problem as long as the potential 
energy features a barrier or transition-state that has to 
be crossed [![. This is the classical Kramers problem Q, 
Fig. [Ha), that has served very well in many areas of 
science. However, the escape-rate problem is known to 
be ill-defined when the particle is trapped in a potential 
well which is the only point of minimum in the potential 
profile (which then either diverges or reaches asymptotes 
along the coordinated axis) The latter case is typi- 
cally exemplified by two noble gas atoms bonded by van 
der Waals forces such that each atom is trapped in the 
Lennard-Jones potential well. The same applies to many 
other sorts of bound states (ej*. diatomic molecules [3,3 , 
deuterons Cooper pairs [J, nuclear neutrons etc.) 
in the absence of long-range repulsion, Fig. [TJd)-(f) . In 
all such cases, the Kramers and transition-state theories 
cannot be applied because the flux driving the particle 
out of the well cannot be defined, and in fact not even an 
energy barrier to be crossed can be identified. This is in 
contrast to many other situations where the competing 
long-range repulsion gives rise to a well-defined barrier, 
Fig. [IJa)-(c). In spite of this theoretical difference, in all 
these systems the observed time of escape is still finite 
and a theoretical estimate is desirable for many appli- 
cations. This is a long-standing and well-known prob- 
lem The common remedy to this difficulty, so far, was 
to employ fictitious absorbing boundaries or long-range 
repulsions to artificially create a barrier or transition- 
state which would allow one to apply the Kramers theory 
for the escape rate over such a barrier Q. This proce- 
dure has the shortcoming that the resulting barrier or 
transition-state is completely arbitrary, and so are the 
results for the calculated dissociation rates. 

Inspired by Peierls 0], here we propose a solution to 
this problem by considering the role of entropy in the es- 
cape process. It is known that entropic contributions can 
influence kinetics in several contexts E3. As an exam- 
ple we recall the problem of the equilibrium flux through 
bottlenecks with a cross-section fluctuating in time stud- 
ied by Zwanzig [ll| , where even in the absence of energy 
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FIG. 1: Examples of situations where the dissociation prob- 
lem is well defined (a-c) and where it is not (d-f). (a) 
The Kramers escape problem Q]. (b) Field-emission of elec- 
trons H, where the external electric field gives rise to a 
well-defined barrier through which the electron can escape 
or tunnel away, (c) The escape of a-particles from atomic nu- 
clei 0], where the repulsion is between positive charges, (d) 
The internuclear potential of a diatomic molecule or molec- 
ular complex [i| . (e) Potential of electrons at metal surfaces 
(thermionic effect) or by neutrons evaporating from hot nu- 
clei [f|. (f) The attractive potential between the two electrons 
of a Cooper pair in real-space [f| . 



barriers due to conservative potentials, the passage time 
is controlled by entropic barriers due to the fluctuations. 
However, the question about the fundamental mechanism 
by which entropy comes into play has remained largely 
unanswered and an analytical framework of general ap- 
plicability is lacking. Here we show how entropy drives 
the flux responsible for dissociating the bound states in 
the absence of stabilizing barriers and analytically calcu- 
late the microscopic dissociation rate of classical bound 
states. The generic theory is validated by showing that 
it is able to recover the classical Saha ionization degree 
in the dilute limit. Further, our approach is fully micro- 
scopic, which allows us time to include the density effects 
in the recombination process. 

Effective pair-interaction energy for dissociation. We 
start with two thermal particles bonded by an attractive 
interaction where it is evident that upon moving apart 
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from one another along the radial coordinate the two 
particles will be less favored in terms of potential en- 
ergy but they will gain a larger free volume, and hence 
more entropy. Furthermore, one should remark that the 
"pair-interaction energy w(r)" between two molecules 
or particles separated by a distance r is usually iden- 
tified by the force / acting between the two particles via 
/ = — ditf(r)/dr; hence one finds the work that can be 
done by the force. This means that w(r) is actually the 
free energy or the available energy for the two-particle 
system J9| because the local collisional physics allows 
us to define an entropy in addition to the potential en- 
ergy. Rigorously, the effective pair-potential for the two- 
particle system is thus given by: w(r) = U{r) — TAS{r), 
where U is the conservative potential energy of interac- 
tion between the two particles and AS is the mixing en- 
tropy of the system. Let us interpret the dissociation co- 
ordinate r as a time-averaged position where the average 
is taken over a time t > t c where r c is the collision time- 
scale between the two bound particles. Then, the contact 
force f c due to collisions between the two particles con- 
strained to remain at close contact over a finite amount of 
time obeys the scaling relation [l4|: f c ~ fcgT/r. For two 
hard-spheres constrained into a small portion of space 
where they collide repeatedly with each other, the con- 
tact force can be integrated to give an associated pair- 
interaction free energy which is related to the entropy of 
the two-particle system: TAS(r) ~ /csTln(r/i?) where 
we chose the integration constant equal to — lni?. In the 
presence of a conservative interaction U(r) the total in- 
teraction force is thus given by / = f c — (dll(r)/dr) and 
the corresponding total effective pair-interaction is given 
by w(r) ~ U(r) — ksT\n{r / R). A more precise form for 
this effective pair interaction, including prefactors, can 
be derived using a different method, which makes use 
of the Onsager excluded volume theory With this 

method one obtains: 
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where V x {r) = 4tit 3 /3, u = 4tt(2.R) 3 /3, and Up is the 
volume of one particle. This equation expresses that fact 
that upon moving the particles apart along the outward 
radial coordinate r there is a net entropy gain which 
arises from the increased number of degrees of freedom 
explored by the particles. Hence the interaction energy 
can be rewritten as 



w(r) = U(r) - k B T 6 In — + 4 - 2 In 2 
V it 
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Based on these considerations, the effective pair- 
interaction potential given by Eq. ([2]) is coarse-grained in 
the sense that is valid only for particles that have been 
kept at close distance for a time t 3> r c [lij . As such, 
Eq. ([2]) can be applied to the dissociation of pairs of parti- 
cles which are bonded by some attractive potential U(r). 
On the other hand, for two particles which approach each 
other from far apart, their mutual interaction prior to 
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FIG. 2: Schematic of the rectangular-well potential around a 
particle of radius R used for U(r) in the derivation, shown by 
the dashed line. It is the effective potential (free energy) w(r), 
solid line, with its entropic correction making the higher sep- 
aration r favorable, that affects the dissociation of the bound 
state. 



colliding happens on a time scale t < r c , i.e. they can- 
not explore their mutual excluded volume. Therefore for 
the recombination process there is no entropic effect and 
w(r) — > U(r). This can be understood by recalling that 
entropy in classical systems is ultimately related to colli- 
sions which are important in the bound-state where the 
particles preserve "memory" of each other [l5{ over the 
life-time of the bound state and this makes their con- 
figurational entropy depend on their separation. In the 
recombination process, however, since the conservative 
potentials U that we consider here do not have a bar- 
rier, the recombination of diffusive particles is diffusion- 
limited and occurs at the very first collision event be- 
tween the two particles. Therefore, no entropic contri- 
bution applies to recombination processes of this kind 
where the recombination rate is uniquely determined by 
the diffusion process in the field of the conservative po- 
tential U. In the following we will use these arguments to 
derive microscopic dissociation and recombination rates 
with the tools of statistical mechanics. 

Dissociation rate. Without loss of generality, but to 
simplify the algebra and emphasize the qualitative point , 
let U (r) in Eq. © be a rectangular- well attractive poten- 
tial, of range 5 (Fig. [2j dashed line) , so that two particles 
constrained within a distance r = 2R + S. Due to the 
entropy contribution, w(r) from Eq. © decreases loga- 
rithmically at large distances, while having little effect on 
the proximity well. Such an effective potential (free en- 
ergy) features a now well-defined barrier over which the 
bound particle can escape. Using the Kramers escape 
theory for the 3D spherically-symmetric problem [l[ , the 
escape time, i.e. the reciprocal of the escape rate re, is 
given by 
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The upper limit C 3> R + 8 is some arbitrary point far 
away along the radial axis and r m j n = 2R + 5 is the min- 
imum of the effective potential w(r); D is the diffusion 
coefficient for this problem. Using Eq. © in Eq. Q the 
integrals can be evaluated analytically which gives the 
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compact expression for the dissociation rate: 
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The escape rate is directly proportional to the diffusion 
coefficient, as in the classical Kramers theory. It is im- 
portant that this dissociation rate decreases with the in- 
creasing width of the attractive well, 5: this is a natural 
reflection of the fact that the effective frequency of the 
particle in the potential well (related to the 'rate of at- 
tempts' to overcome the barrier) increases as S — > 0. The 
height of the effective energy barrier in the exponential 
thermal-activation term, is given by A/k B T, that is, by 
the depth of the original square well potential. This ex- 
plains why a number of ad hoc theories produce results 
that are qualitatively valid; the crux here is in the detail. 

Recombination rate. Once the dissociation of the 
bound state occurs, i.e. the bound particle has crossed 
the barrier in the effective potential of Fig. [2 the two 
particles move apart. For a dilute thermal system, the 
mechanism of recombination back into the bound state 
is controlled by the diffusive transport. With interacting 
particles, the recombination rate is determined by solving 
the stationary Smoluchowski (diffusion) equation in the 
field of force due to the interaction. For the attractive 
square- well case the solution is: 
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Since 5 is finite and for A/k B T > 1, the first term in the 
denominator is small compared to the second term due 
to the exponential factor, and thus we obtain: 
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This result can be extended to denser fluids using 
a generalized diffusion coefficient given by D{n) — 
(D/k B T)dH(n)/dn where II is the osmotic pressure of 
the fluid medium and n = N/V is the average number 
density of diffusing particles [13]. Using this D(n) in 
the above solution to the Smoluchowski diffusion equa- 
tion, one obtains the generalization of the recombina- 
tion rate for dense and crowded systems [f7j : n rec 
4irD(2R + 5)H(n) / nk B T . The correction can be easily 
evaluated once the equation of state (EOS) of the fluid 
is known. In the simplest non-ideal case one can express 
the EOS in terms of the virial expansion [l8( . For model 
hard-sphere fluids one can use the Carnahan-Starling 
EOS 0: U/nk B T = Z($) = (1 + + (f) 2 - </> 3 )/(I - 0) 3 , 
where <f> is the fraction of occupied volume in the system, 
e.g. 4> — (47r/3)i? 3 n. This leads to the recombination- 
rate coefficient K rec = 47tD(2R+S)Z(4>), which would be 
relevant for a crowded environment. For a dilute solu- 
tion of biomolecules which bind to form complexes (e.g. 
receptor- ligand binding) in a sea of other particles of sim- 
ilar size with which they interact by steric repulsion, this 
leads to recombination rates several times larger than the 



Smoluchowski estimate. We remark that for ^ > 0.40, it 
is no longer safe to neglect entropic effects on recombi- 
nation, due to the short-range liquid-like structuring (in- 
homogeneous local density) around the reactants which 
cause the incoming particle to collide with an increasing 
amount of crowders before colliding with the second re- 
actant. This clearly leads to an entropic repulsion upon 
recombination, as discussed in [l7| . These considerations 
fix the upper validity limit of our theory to <f> ~ 0.45. 
Extension to glassy systems is non-trivial and will be the 
object of future work. 

If the dense fluid is made of charged particles, such 
as in high-density plasmas, then the transport is no 
longer diffusive but drift-controlled. In this case the 
oppositely charged particles attract each other via the 
Coulomb force ~ (Ze) 2 /r 2 , where Ze is the electric 
charge of the particles. The crossover from diffusion- 
controlled to drift-controlled transport typically happens 
when n[(Ze) 2 /k B T]a ^> 1, where a is the collision cross- 
section. In this case the total current of oppositely- 
charged particles entering the interaction volume is given 
by: 4irr 2 Vdn, where Vd = (Ze/r 2 )(/i_ + /i + ) is the drift 
velocity in the Coulomb attraction field of approaching 
particles and (i the electrical mobility. The electron-ion 
recombination-coefficient for dense plasmas is thus given 
by Krec ~ 47re/x_, assuming that the electron mobility 
fi_ is much greater than the ion mobility and Z = 1. 
This is a well-known result from the kinetic theory of 
plasmas [20l |. 

Dissociation equilibrium. The equilibrium of asso- 
ciation and dissociation reactions is of extreme impor- 
tance for in vivo physiological processes, e.g. enzy- 
matic activity based on receptor-ligand complexes or 
protein aggregation [111, [22| . We will give a treatment 
of the general case of two classical particles bound by 
an attractive potential which can be directly applied 
to biomolecules. The net attractive potential U(r) is 
of the kind of Fig. [IJd). Given these conditions, the 
chemical equilibrium constant K for the dissociation- 
recombination process B + C ^ BC is given by the law 
of mass action as 
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If N = N B c + N B = N BC + N c is the total number of 
each species present in the system, the total number of 
particles, both bound and dissociated, is given by N — 
N (l + a), where a = N B /N a = N c /N is the "degree of 
dissociation" [24]|. The relations for the molar fractions 
of the various species follow at once: N B /Nq = a/(l + a) 
and N BC /N Q = (I - a)/(l + a). 

Using Eq. (f4]) and the high-density extended Smolu- 
chowski rate for Kdiss and n rec , respectively, this leads 
to the following expression for the dissociation degree in 
equilibrium: 
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In the limit of long-range attraction, 8 3> 2R, this sim- 
plifies further and we get the following expression for the 
degree of dissociation a: 
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In biological systems one often encounters situations 
where the interaction range can be significant, such as 
in the ubiquitous case of hydrophobic attraction p5| . 
Equation ([9]) shows that the dissociation degree in cases 
of biological relevance can have a sensitive dependence 
upon the interaction range 8, which has been neglected 
in previous theories. Furthermore, Eq. ([9]) takes into ac- 
count the enhanced recombination in crowded systems. 
In the opposite limit of short-range attraction (sticky or 
adhesive particles [26|), 8 <C 2R, one obtains: 
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with a much weaker dependence on 8. Finally, we 
should remark that in the above treatment we assumed 
n B = nc, which describes dissociation and ionization ki- 
netics. In receptor-ligand kinetics, however, the concen- 
tration of ligands normally overwhelms the concentration 
of receptors, i.e. nc S> n B . Since we derived intrin- 
sic rates, we can apply our theory to any kind of kinet- 
ics, including the pseudo-first-order situation nc n B . 
In this limit, the kinetics is of the Langmuir type (23| : 
9 = Kn c /(l + Kn c ), where 9 = N BC /(N B + N BC ), and 
K is still given by Eq. (J7J with the rates for dissociation 
and association derived here. 

Ionization equilibrium in plasmas. We can further 
test this theory by addressing the problem of the ther- 
mal dissociation of atoms (ionization). This problem was 
famously addressed in the 1920's by Saha [24| by com- 
bining the statistical mechanics of the ideal gas with the 
chemical equilibrium (detailed balance) assumption. By 
considering the chemical equilibrium of the ionization re- 
action A ^ I + + e~ , the chemical equilibrium constant of 
the reaction is given by the same expression as Eq. (0. 
In the dilute limit, Eq. @ for K rec applies, and using 
Eq. ([¥]) for Kdiss we obtain 
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where N = Na + Nj+ + N e ~ . The interaction volume, 
as we have shown above, is given by v = [{2R + 8) 9 — 
(2R) 9 ]/(2R + 8) e . The extension of the radial probabil- 
ity amplitude for a bound electron, A e - , is approximately 
given by the width of the attractive well, and therefore 
we take 8 ~ A e - . In the classical approximation valid at 
high temperatures, A e - is the thermal de Broglie wave- 
length, A e - = %/ yJm e k B T '/2tt. Furthermore, one has 
A e - > A 7 +, such that v = [(A/+ + A e - ) 9 - (A/+ ) 9 ]/(A /+ + 
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i.e. the same limit as used in Eq. ([5]). 



With this replacement in Eq. our theory for a di- 

lute plasma where free diffusion is the main transport 



mechanism correctly reduces to the Saha equation [24|, 
apart from a numerical pre-factor of order unity. This re- 
sult demonstrates the validity of the proposed calculation 
scheme. 

The Saha equation works well for dilute plasmas, but 
it breaks down in the core of the stars [27|. In fact, the 
opacity in the interior of stars (such as the Sun) is ex- 
tremely low for that high density, and this implies that 
the interior is formed by completely ionized matter with 
a ~ 1. Evaluating the Saha equation at the typical con- 
ditions of the stellar interiors gives a ~ 0.7 — 0.8 that 
is incompatible with the observed vanishing opacity. We 
can now extend the theory to the high-density regime 
because we can include transport effects, that are impor- 
tant at high density, via the recombination rate which 
enters in our model. Using the recombination rate for 
high-density plasmas given by K rec « iirefi- in the equi- 
librium constant for atomic dissociation K — Kdiss/K rec , 
we obtain 
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where we used the Einstein relation eD/k B T = fi for the 
single-electron ionization, Z = 1, and continued using 
8 ~ A e - for the bound state of electron. Let us calcu- 
late a for the conditions corresponding to the core of the 
Sun, i.e. T = 10 7 K, N/V = 10 26 cm~ 3 , and assume pure 
hydrogen. The Saha equation gives a ~ 0.79, whereas 
from Eq. (fT2j) we obtain a much more realistic a ~ 0.99. 
Clearly, Eq. (fT2j) is in excellent agreement with the exper- 
imental data indicating fully- ionized plasma conditions in 
the core of the Sun [28| • 

To conclude, we have proposed a theoretical scheme 
that allows one to calculate the intrinsic dissociation 
rate of bound states with purely attractive potentials, 
something so far possible only for potentials that feature 
long-range repulsion competing with attraction. In all 
the other cases, including dissociation of neutral atoms, 
molecules, Cooper pairs, molecular complexes, this was 
not possible and required various ad hoc assumptions. 
We validated our theory for its predictions of the ioniza- 
tion degree of hydrogen-like atoms, a problem relevant 
in several branches of astrophysics. Our theory correctly 
reproduces the Saha formula at steady-state in the ideal- 
gas limit. Further, it can be easily extended to deal with 
the high-density limits. In fact it predicts full ionization 
in the core of stars, thus dramatically improving over the 
Saha equation that predicts partial ionization. Further 
applications are in biology where it allows one to calcu- 
late dissociation equilibria of receptor-ligand or protein 
complexes in presence of cellular crowding. 
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